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The use of electric propulsion is more prevalent than ever, with industry pursuing all electric 
orbit transfers. Electric propulsion provides high mass utilization through efficient propellant 
transfer. However, the transfer times become detrimental as the AV transitions from near- 
impulsive to low-thrust. Increasing power and therefore thrust has diminishing returns as the 
increasing mass of the power system limits the potential acceleration of the spacecraft. By 
using space-to-space power beaming, the power system can be decoupled from the spacecraft 
and allow significantly higher spacecraft alpha (W/kg) and therefore enable significantly 
higher accelerations while maintaining high performance. This project assesses the efficacy 
of space-to-space power beaming to enable rapid orbit transfer while maintaining high mass 
utilization. Concept assessment requires integrated techniques for low-thrust orbit transfer 
steering laws, efficient large-scale rectenna systems, and satellite constellation configuration 
optimization. This project includes the development of an integrated tool with 
implementation of IPOPT, Q-Law, and power-beaming models. The results highlight the 
viability of the concept, limits and paths to infusion, and comparison to state-of-the-art 
capabilities. The results indicate the viability of power beaming for what may be the only 
approach for achieving the desired transit times with high specific impulse. 
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Nomenclature 

Area, m 2 
Semimajor axis 
Eccentricity 

Acceleration due to thrust (m/s 2 ) 

Specific angular orbital momentum (m 2 /s) 

Inclination 
Semilatus rectum (m) 

Radial distance from central body (m) 

Proximity quotient (s 2 ) 

Scaling function 

Weight for orbit element oe 

Orbit element (one of a, e, i, co, O) 

Angle of the projection of the thrust onto the orbit plane with respect to the local horizon, positive 
outwards (rad) 

Angle between the thrust and the orbit plane, positive in the direction of the oribital angular 
momentum (rad) 

True anomaly (rad) 

Gravitational parameter of the two-body system (m 3 /s 2 ) 

Argument of periapsis (rad) 

Longitude of the ascending node (rad) 

Maximum over true anomaly of the maximum over thrust angle, or giving said maximum 

Component along the orbital angular momentum vector 

Radially outward component 

Circumferential component 

Target value 


* Mission Analyst, Technology Development and Transfer Office, ZP30, NASA MSFC. 
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I. Introduction 


W hile electric propulsion is enabling industry to maintain cost competiveness, progress is only evolutionary and 
is limited by the power system mass, efficiency, and cost. Solar cell efficiencies continue to increase at a 
reasonable cadence. New solar array packaging and deployment methods should also yield continued evolutionary 
progress in specific power in the near-term. Even with the trends of progress, they are only evolutionary and do not 
alone solve the limiting detriment to electric propulsion, long transit times. 

The transit time is a function of the 
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spacecraft acceleration and therefore specific power; 
power available for propulsion (converted to thrust) 
divided by the mass of the spacecraft. While the 
power system efficiency and absolute power 
capabilities are growing, the desired spacecraft 
masses are also increasing. As shown in figure 1 , the 
trend is favorable to electric propulsion and 
exceeding 4 W/kg, but still necessitates relatively 
long transit times over chemical alternatives. The 
results shown in figure 2 indicate the required 
thrusting times departing from a geo-transfer orbit 
(GTO) and low earth orbit (LEO) at near-term 
specific powers using electric propulsion. 

Figure 3 illustrates the thrusting times 
achievable from GTO-to-GEO and LEO-to-GEO for 
significantly higher specific powers. Results in 
figures 2 and 3 were generated using NASA’s 
SEPSPOT tool. The transfer vehicle would still require payload power, the propulsion system, a power receiver, and 
thermal radiators, but a basic implementation of the power beaming spacecraft should yield an immediate leap in 
capability from 4W/kg to 40W/kg with a modest electric propulsion system. For a 2,000s Hall system, GTO-to-GEO 
transits would go from 118 days of thrusting to only 12 and LEO-to-GEO transits would decrease from 235 days to 
24 days. Mass fractions from GTO and LEO would increase from 57% to 87% and 27% to 75% over chemical 
bipropellant respectively. If successful, the paradigm of launching large payloads to GTO for chemical insertion 
would likely shift to mass efficient lower cost launches to LEO. 
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Figure 1. Trend of increasing COMSAT power and 



Figure 2. Near-term GTO-to-GEO transit times. 


Figure 3. Transit times with SSPB. 


II. Concept Description 

The concept proposed is fundamentally high specific power (W/kg) electric propulsion enabled by disassociating 
the power generation from the transfer vehicle. The VASIMR engine made headlines with claims of orbit transfers 
from Earth to Mars in less than 40 days; truly revolutionary and high value to NASA. However, upon closer look; it 
was not the propulsion system that was revolutionary. Conventional electric propulsion thrusters can achieve 
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comparable performance for the Earth-to-Mars transit 
provided they use the same assumption for power 
system specific power, >2,000W/kg. The specific 
power was revolutionary, although NASA is not 
investing in the required power system to achieve 
those goals. 

One method to achieve extremely high specific 
power is to remove the power production from the 
spacecraft entirely. Power generation spacecraft can 
be placed in strategic orbits and beam the power to 
the transfer vehicle when line-of-site (LOS) exists 
and power generation spacecraft are illuminated. 

This concept is illustrated figure 4. Note that power 
receiving and thermal management systems are still 
required on the transfer vehicle. 

III. Power Transmission Analyses 

m i i Figure 4. Illustration of space based power beaming. 

The emphasis ot the power system analyses is ° 

focused on the transfer vehicle to determine the potential power-to-mass ratios achievable for the transfer vehicle. 
First, power is acquired on the power-beamers through the use of solar arrays. This DC power is then converted to RF 
power so that it may be transmitted to the rectenna (rectifying antenna) on the transfer vehicle. A rectenna captures 
incident RF power and transforms it to DC power by a diode based converter. An HG filter ensures impedance match 
between the rectifier and the antenna for optimal power transfer. The output DC filter smoothens the output DC voltage 
and current by attenuating high frequency harmonics present in the RF signal or generated by the highly nonlinear 
rectification process 1 . 

A. Power Conversion 

The most limiting factor of the power-beaming concept is the efficiency of the power conversion process. Table 1 
shows the efficiencies used for this study, based on previous GEO-to-ground wireless power transmission studies 1,4 . 
The DC-RF converter efficiency was determined by assuming the use of a Gyrotron. Gyrotrons enable high power 
transmission at 94 GHz with efficiencies up to 60%. Collection efficiency was based on the divergence of the beam 
over distance 2 . Due to the large divergence of the beam, the efficiency is poor at long distances since the receiver is 
not large enough to capture much of the energy. Rectenna efficiencies were determined based on the state of the art 
efficiencies for rectenna arrays. The conversion efficiencies were scaled for 94 GHz by calculating the effective area. 
As the frequency increases and the effective area of the rectenna decreases, the power density required to achieve a 
given efficiency increases. An efficiency curve was created based on McSpadden’s and Mankins’ data 3 for 2.45 GHz 
and then scaled to the appropriate levels for a 94 GHz transmission. This gave the conversion efficiency based on the 
power density. The power density was calculated by using a series of receiver sizes to determine the power density 
in concentric rings within the antenna. 

Though this process may seem fairly inefficient, a rectenna yields a specific power of 4000 W/kg 3 , while the use 
of on-board solar arrays would yield a specific power of approximately 70-90 W/ kg 4 . 


Table 1. Summary of constant efficiencies used for study. 


Sources of Inefficiency 

Efficiency 

Notes 

RF Circuit Efficiency 

DC-RF converter efficiency 

0.600 

Assumes gyrotron converter efficiency 

RF filter insertion loss (IL) 

0.891 

Estimated total IL = 0.5 dB 

Transmitting Antenna Efficiency 

Subarray random electronic failures 

0.960 

Estimated 2% failures 

Meteorite hit element failures 

1.000 

100 failures/year 

Amplitude error 

0.996 

±1 dB amplitude deviation 

Phase error 

0.978 

±15° phase deviation 

Antenna aperture efficiency 

0.980 

Conductive losses in aperture 

Transmitter scan loss 

1.000 

Assumes broadside radiation 

Mismatch loss 

1.000 

Assumes Array VSWR =1.0 
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Beam Coupling Efficiency 

Collection efficiency 

Varies 

Depends on distance between transmitter and receiver 

Polarization loss 

1.000 

Assumes near perfect alignment 

Rectenna Efficiency 

Rectenna random failures 

0.990 

Estimated 1% failures 

RF filter insertion loss (IL) 

0.891 

Estimated IL = 0.5 dB 

Rectenna scan loss 

1.000 

Rectenna tilted to avoid loss 

Mismatch loss 

1.000 

Assumes Array VSWR =1.0 

Rectenna efficiency 

Varies 

Depends on distance between transmitter and receiver 

DC-DC converter efficiency 

0.950 

Assumed dc-dc converter efficiency 

System efficiency 

<41.05% 

Not including collection and rectenna efficiencies 


B. System Optimization 

Figure 5 shows that as the 
frequency of the wireless power 
transmission increases, the divergence 
of the beam decreases. This allows for 
improved efficiency of the power 
collection with the same size antenna, 
or alternatively, the same efficiency 
while using a smaller receiver or 
transmitter. In this case, the size of the 
receiver limits the utility of this 
technology. Therefore that is where it 
would be most advantageous to reduce 
the size of the receiver. Increasing the 
frequency also improves the 
conversion of RF power back to DC 
power, since the rectenna is more 
efficient when operating at higher 
power densities. 

However, use of higher frequencies 
is not without its downsides. One limitation to higher frequencies is the decreasing efficiency of solid-state DC-to- 
RF converters at higher frequencies. For example, at 4 GHz, the efficiency of a converter is 85%. However, at 30 
GHz that efficiency falls to only 65%. Where they exist, alternative technologies which provide higher efficiencies 
at higher frequencies are currently not space ready. Another very important limitation to higher frequency 
transmission is the ability to build large receivers which are smooth enough to properly receive the signals. The 
smoothness of the receiver surface needs to be on the order of the wavelength being received. Therefore at higher 
frequencies, the precision needed in the surface of the receiver increases and this is limited by the size of the receiver 
technologies. 

C. Transmission over Large Distances. 

Over short distances, the efficiency of the rectenna and the collector can approach 100%. Due to system losses 
and DC-to-RF efficiency, the total transmission efficiency is nearly 41%. Although the rectenna conversion and 
collector efficiencies decrease significantly at longer distances, the distance losses can be greatly reduced by placing 
the power generation satellites at strategic locations. For example, in order to raise a GTO orbit to GEO, the perigee 
needs to be raised which requires maneuvers at apogee. Therefore, placing a power generating satellite at GEO at 14° 
inclination will bring the satellites as close as 10,500 km throughout the orbit transfer which will give a power 
collection efficiency of 84%. Therefore, the highest thrust can be achieved at optimal locations in the orbit leading to 
faster orbit maneuvers with only a few satellites in key locations. 

IV. Performance Assessment Methodology 

There are several methods available for calculating optimal or near-optimal low thrust trajectories. The options 
considered for this project were Q-Law, Kluever’s method, Betts’ method, and FASTSTOP. Q-Law was ultimately 
chosen based on its ability to target all five orbital elements, allowing for use with a wider range of orbital maneuvers, 


Aperture Equation Results 



Figure 5. Aperature Equation Results; as the frequency of the wireless 
power transmission increases, the divergence of the beam decreases. 
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and the ability to know that satellite’s position exactly throughout the trajectory which aided in calculating the power 
received at the satellite. While the method used by Betts could provide an optimal solution to the problem, the 
optimization method is more complex to implement and should yield results within a few percent of the Q-law 
approach; sufficient for a preliminary concept viability study. The other two options, Kluever and FASTSTOP both 
use orbital averaging which results in faster trajectory simulations, but results in a loss of exact position information 
which is desired to assess the dependence of the orbital location of the power transmission spacecraft. 

A. Q-law 

Q-law is a Lyapunov feedback control law which determines the best thrust direction at the current point in the 
orbit based on a promity quotient Q which as defined as: 

Q = (1 + WpP) Yoe W oe S oe [ d(0 !’° er) | 2 , for oe = a, e, i, m, H (1) 

L oe XX J 

where the classical orbital elements of semimajor axis, eccentricity, inclination, argument of periapsis, and longitude 
of the ascending node are used to target an orbit but not an exact position. This expression gives a measurement of the 
distance from the target orbit. The W ’ s are weighting values, the S’s are scaling functions, and d(oe, oe T ) is a distance 
function indicating how far the current orbital element value is from the target value. P can be used to apply a penalty 
function while oe xx are analytical expressions for the maximum rate of change of the orbital elements in a given orbit. 

Taking the derivative with respect to time of the quotient, Q, gives the rate at which the orbit is approach the target 
orbit. This value is given by: 


dQ 

dt 


dQ 
' doe 


oe 


( 2 ) 


It is a function of the angle at which the thrust is applied since the oe are Gauss’s form of the variational equations 
for the orbit and those equations are functions of the thrust angles. Therefore, to find a near optimal low thrust orbit 
transfer, the time rate of change of Q is minimized all along the trajectory. 

The angles a and (3 produced during the orbit simulations matched the form of the thrust angle evolution as 
derived by Anastassios Petropoulos 5 for the case of a GTO-to-GEO or LEO-to-GEO orbit transfer. Similar 
performance was also seen with most over orbit maneuvers. However, the implementation of Q-law was unable to 
successfully perform transfers which involved more than a 90° change in inclination and could not be compared 
against the GTO-to-Molniya case presented by Petropoulos. The implementation of Q-law was also compared 
against the optimal time transfer produced by the low thrust mission design tool SEPSPOT, allowing for a rough 
check of the trajectory times which were approximately 10% longer with Q-law than with SEPSPOT. 


B. IPOPT 

IPOPT is a free optimization tool with a MATLAB interface which was used in combination with Q-Law to 
determine the low-thrust trajectories. Specifically, IPOPT was used to determine which angles minimized the 
derivative of Q at a specific point in the orbit. Those angles are then chosen as the thrust angles until the next 
optimization is performed. 

IPOPT uses an interior point line search filter method to find a solution to the optimization problem. This method 
requires not only a calculation for the objective function, in this case Q ', which is being minimized but also the gradient 
and the Hessian matrix for the objective function with respect to the variables being optimized. The form is the 
Hessian matrix associated with this optimization problem is shown in Equation 1 . If constraints are used, then those 
constraints and the associated Jacobian matrix also need to be provided. In the case of the angle optimization for Q- 
Law, the angles had bounds but no constraints so no Jacobian matrix was required. 


' d 2 Q 

d 2 Q ' 

da 2 

dad/3 

d 2 Q 

d 2 Q 

Ldpda 

Jp 2 . 
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V. Results 


Three missions were evaluated, including a GTO-to-GEO and LEO-to-GEO transfer assuming the use of the GPS 
constellation as power beamers and a LEO-to-GEO transfer again using the GPS constellation, but with an additional 
four power beamers placed in LEO. 

C. GTO-to-GEO 

The first case considered was a 
transfer orbit from GTO-to-GEO 
where the initial orbit had an apogee at 
42,164 km and a perigee at 6,521 km 
with an inclination of 28.5°. For this 
orbit maneuver, the transfer times 
ranged from 7 days to over 700 days 
depending on the mass of the satellite, 
the power being transmitted and the 
maximum allowed received power. 

The results for 50 kW transmitters on 
the power generation satellites are 
shown in Figure 6. As expected, as 
satellite mass increases, transfer time 
increases. Furthermore, combinations 
of higher power transmitters and 
receivers both contributed to lower 
transfer times, however, the best 
transfer times were achieved when the 
maximum power on the transmitter 
and receiver were close in value. This 
can additionally be seen in Figure 6, which shows transfer times for various combinations of transmitted and maximum 
allowed received power all assuming a satellite mass of 5,000 kg. Figure 7 also shows that when the transmitted 
power is less than the receiver’s maximum, increasing the receiver power capibilites will have no impact on transfer 
time. Similarly, when receiver’s maximum allowed power is less than the transmitted power, increasing the 
transmitted power will have little effect on transfer time since the transfer vehicle is already receiving its maximum 
power. 

D. LEO-to-GEO 

The second case considered a transfer 
orbit starting in LEO, specifically a 
circular orbit with a semimajor axis 6,671 
km and an inclination of 28.5°. Since this 
starting orbit is at a lower altitude, the 
transfer requires a larger delta-V and, as 
expected, takes longer to complete. Figure 
4 shows that the transfer times for the 
LEO-to-GEO cases ranged from 9 to 900 
days, taking about 20% longer than the 
GTO-to-GEO cases. 


Summary Point Cases: GTO-to-GEO 
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Figure 7. Summary point cases for a transfer between GTO and 
GEO at a satellite mass of 5,000 kg comparing different transmitted 
powers. 



Figure 6. Transfer time versus satellite mass for 50 kW transmitter 
maximum power for select receiver maximum power for GTO-to-GEO. 
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E. LEO-to-GEO with Additional 
Four Beamers 

The last case that was tested was 
the same transfer orbit as the LEO-to- 
GEO transfer in Case B above, 
however, an additional for beamers 
were placed in LEO. Figure 6 shows 
the transfer times for this case. For 
both maximum allowed receiver and 
transmitted powers of 50 kW, the 
LEO-to-GEO case with only GPS 
beamers yielded transfer times 
ranging from 18 to 185 days, where 
the case with the additional four 
beamers gave transfer times ranging 
from 18 to 178 days. The addition of 
the beamers enabled the transfer 
vehicle to have access to a closer 
power source at the start of its 
transfer, thus decreasing the transfer 
time from the GPS-beamer-only case 
above by about 4%. 



Figure 8. Transfer time versus satellite mass for a LEO to GEO transfer for 
50 kW transmitter maximum power for select receiver maximum power. 



Figure 9. Transfer time versus satellite mass for 50 kW 
transmitter maximum power for select receiver maximum 
power. 


Caie B & Case C Comparison 



Figure 10. Transfer time versus satellite mass for two power 
beamer constellations both with 50 kW transmitter and 
receiver maximum power. 


VI. Work to Go 

Although tests have been run using the GPS constellation as the power generation satellites, tests will also need to 
be performed using considering other cases such as power beaming satellites in LEO or GEO orbits. Also, although 
the tool currently has the capability to perform brute force optimizations for the position of consetllations of satellites, 
it would also be useful to implement a true optimization of the satellite positions. With the current tool, a constellation 
optimization was attempted, however, since optimizing the constellation requires calculating the transfer time with 
each position for the power generating satellites, the optimization for a GTO-to-GEO orbit transfer took 4 to 8 hours 
and produced values which appeared to be local, not global, minima. If the constellation optimization be performed 
in a timely manner, then the transfer times could likely be reduced by 10-20% if not more. 

During the implementation of the power transmission between satellites, the interference between power beams 
coming from different satellites was not fully investigated. The power received was was assumed to add linearly. 
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This is indeed the case when the beams are not phase synchronized. However, when receiving power from multiple 
sources if the phase of the transmitters can be synchronized, then the power received from N satellites can be up to N 2 
the power received from one satellite. Although the best case scenario occurs when the beams are exactly 
synchronized, most of the extra power can be attained as long as the beams are synchronized to within one-eighth of 
a wavelength. It may be useful to consider this effect in future space based power transmission work. 

Additionally, the rectenna conversion efficiencies used with the calculation of power transmission were based on 
covering an entire antenna surface with the rectenna array. Considering that rectenna efficiency increases with higher 
power concentrations, an alternative method might be to concentrate the signal into a smaller rectenna area. This 
could be done by using the receiver as a reflector which sent the signal to a smaller rectenna array, much like radio 
telescopes use dishes to concentrate the signal. This method could be used to both increase the overall efficiency and 
to potentially reduce the required mass of the receiver 

VII. Comments on Architecture Implementation 

This preliminary study assessment the potential to decouple the power generation system from orbit transfer 
vehicles to enable a dramatic improvement is spacecraft power to mass ratios beyond the current evolutionary path. 
One of the challenges with revolutionary concept implementation is the cost to realize return on investment. 
Momentum exchange Electrodynamic Reboost (MXER) and space elevator tethers are concepts to change space 
transportation architectures, but require tremendous capital investments prior to first application. The use of space- 
to-space power beaming, however, offers a practical infusion path. The system can be analogous to the US investment 
in roads. The US did not build a highway across the country as the first road. Instead, local roads were established, 
then intrastate and eventually interstate highway systems. It would not have been practical to begin with billions in 
investments in roads as an initial investment, however, the long duration investments have led to a transportation 
infrastructure critical to logistics of the nation. Space-to-space power beaming can similar be implementing in 
stepping stones. CubeSats are capable of deploying 100m 2 surfaces and a LEO power beaming system could enable 
CubeSat power and electric propulsion for orbit transfer and station keeping without the need for billions of spacecraft 
investments prior to first application. This would allow both a low cost implementation approach while closing a 
critical gap of SmallSat technology. Adding power beaming functionality to GEO spacecraft can be used to support 
GTO to GEO orbit transfer and over time, a power beaming infrastructure could evolve to as a critical asset to the 
logistics of geocentric space transportation. 


VIII. Summary 

The implementation of a geocentric power beaming architecture can revolutionize the fundamental methodology 
of placement and transfer of space assets. Leveraging advanced in-space propulsion is known to reduce mass required 
to LEO and therefore cost of access to space. Electric propulsion systems have always been power limited, and the 
dissociation of the power production from the transfer spacecraft can create a great leap in capability for reducing 
orbit transfer time; abating its limiting attribute. The power beaming architecture and infrastructure would result in a 
beneficial change to NASA’s long-range plans and benefit nearly all space enterprise. The implementation of this 
concept will yield a revolutionary improvement of performance, mass reduction, and cost reduction of future 
spacecraft. The architecture can also enable SmallSat propulsion to augment the already expanding market by 
increasing capabilities and addressing the need for deorbiting spacecraft with otherwise limited onboard propulsion. 
This concept can leverage NASA and DoD investments in propulsion systems, large deployable arrays, advanced 
thermal management technologies and transmitter and receiver technologies. If fully implemented, the concept 
appears to be the only viable method to achieve both high mass utilization and rapid transfers from launch vehicle 
insertion to final operational orbits. While the concept is new relative to previous power beaming studies, it warrants 
further study to explore the feasibility and viability to meet entrance criteria of future investment opportunities. 
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Appendix A: Steering Law Implementation 

The control law implemented is based on Q-law. Q, which is defined as follows: 

Q = (1 ~ W P P) £ oe W oe S oe [ d(o !' oe ^ | 2 f for oe = a, e, i,co,£l (4) 

L oe XX J 

where classical orbital elements are used and an orbit, not a position is targeted, so true anomaly is not included. W P 
is the weighting for the penalty function which is generally either 0 or 1 . P is a penalty function which causes the 
control law to avoid certain regions, for example in this case as in Petropoulos’ paper, a minimum-periapsis penalty 
function is used. It is defined by the following: 


3 \ r pminJ 


(5) 


P = e 

and in our case k = 100 and r pmin = 6578 km. W oe is the weight placed on driving each orbital element to its 
target value. For this model, all the W oe were set to 1 until the final orbit was reached, except in the case with no 
eccentricity and no inclination, when the argument of periapsis, co, and the right ascension of the ascending node, H 
are not targeted and set to 0. S oe is a scaling function used to prevent non-convergence. It only applies to the 
semimajor axis. Therefore, S oe = 1 for oe = e, i, co, H. However, for the semimajor axis it is: 


1 + 


i 

412 


{ q-fl'A 4 
V 3 a T ) 


Finally, in the equation for Q, d(oe, oe T ) is a distance function indicating how far the current orbit is from the 
desired orbit. This distance function is defined as: 

( oe - oe T foroe = a,e,i 

{oe, oe T ) - [ cos -i[ CO s(oe - oe T )\ for oe = co,£l 
where the distances for co and H are defined so that the shortest distance to the target will always be used. 


( 6 ) 


(7) 


The terms oe xx indicate maximum value of the time derivative of the orbital element out of all thrust angles and 


all true anomalies. These values can be given analytically as follows: 

& xx — 




^rr 


a 3 ( 1+e) 
ju(l-e) 
2 Pf_ 
h 

Vf 




1 cos 2 a)-e|cos a)|) 
Vf 


COy 


xx hsmi^Jl-e 2 cos 2 co-e|sina>|) 

xi = Vp 2 cos 2 e xx + (p + r xx y sin 2 6 X 


(Or 


= Ojcosil 


( 8 ) 

(9) 

( 10 ) 

( 11 ) 

( 12 ) 

(13) 


There are two values of oe xx for co because it contains components in both the in-plane and out-of-plane directions. 
In oe r 


the in-plane and out-of-plane components are weighted in the following manner: 

^xx ~ (^xxi "f 0.01d> xxo )/(1.01) 
while for all the other orbital elements, oe xx 


(^xxi "f 0.01 C) xxo j 

oe xx . Additionally, the values cos 6 XX and r xx are defined as: 


cos 6 XX = 


2e 3 
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+ ■ 


27 


l-e z 

2e 3 
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+ ■ 


27 


r xx = — 

l+ecos 0 X 


(14) 

(15) 

(16) 


To apply Q-law, the in-plane and out-of-plane angles, a and /? respectively are chosen to minimize the rate of 
change of Q, which is given by: 


dQ _ v dQ . 
dt ~^ oe doe° e 

In this equation the values, oe , are the most straightforward, being given by Gauss’s form of the variational 
equations, which are: 

dO. _ rsin(0+co) f 

~Jh 


dt 

di 


h sin i 
r cos ( 0 + 0 )) 


dt h f h 
ft = ld~ v cos 9 fr + (.P + r ) sin 9 fe] 


r sin(0+a>) cos i 
h sin i 


fh 


(17) 

(18) 

(19) 

( 20 ) 
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de 

dt h 
h 

dt r 2 


da 

dt 


= ir{ e sin ( 0 )/r +7/0) 


{p sin 6 f r + [(p + r) cos 6 + re\f e } 
77 = 77 + ^ [p cos ( 0 ) fr-(p + r) sin 0 f e \ 


( 21 ) 

( 22 ) 

(23) 


where p is the semilatus rectum given by p = a(l — e 2 ),h = yj ga( 1 — e 2 ) is the specific orbital angular 

$ ti— g 2 ) 

momentum, and r = 1H _ ecos g is the radius at the current position. The thrust acceleration components are given in 
the radial, tangential, and normal directions as follows: 

f r =f cos /? sin a (24) 

fe=f cos p cos a (25) 

fh = f sin P (26) 

The partial derivative of Q with respect to the orbital elements are given by the following: 

^ = f(l - WpP ) fe oe 2 W oe d Joe^l rt / J. frKoe.oer) _ iCowr) doe xx ^ _ ^ ^oe^l 

doe l y J V-'oe oe doe oe xx \oe xx doe oe xx doe )\ y doe ^ oe oe oe L o'e xx J j 

(27) 

where the sums are over oe = a, e, i, 00 , (1 and all partial derivatives in one equation are with respect to the same 
orbital element. Most of the partial derivatives in the above equation are trivial, however, the nonzero derivatives, 

doexx 
doe 


, can be written analytically as: 
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In order to apply Q-law, a method to calculate the above derivatives and minimize dQ/dt at each point needs 
was implemented. Since the angles were not being changed very rapidly, updates to the thrust angles were made 
after l/54 th of the current period of the transfer orbit, which amounted to the satellite moving approximately 6.66° 
along its trajectory. This spacing was chosen to perform the assessment with practical speed of the trajectory 
calculations while ensuring that the angles generally converged to the correct values. When increasing the 
semimajor axis of the orbit, the in-plane angle, a , was limited to the -n ton radians range and continuously 
increased within that range. The out-of-plane angle, /?, was periodic within its range of — ^ to ^ radians. Any abrupt 
jumps in these angles indicated that the angle had converged to an incorrect value. 
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